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Abstract 

We use Monte Carlo simulations to obtain an improved lattice measurement of the critical cou- 
pling constant [A//i^]^^^^ for the continuum (1 + l)-dimensional (A/4)(/)^ theory. We find that 
the critical coupling constant depends logarithmically on the lattice coupling, resulting in a con- 
tinuum value of [^//^^jcrit ~ ^0-^i o5' considerable disagreement with the previously reported 
[A//i^]^^^^ = 10.26^Q4. Although this logarithmic behavior was not observed in earlier lattice 
studies, it is consistent with them, and expected analytically. 
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I. INTRODUCTION 



The two-dimensional field theory specified by the Euclidean Lagrangian 

CE = lwr+li^i<p'+^^' (1) 

exhibits a phase transition between a symmetric phase with (6) = and a phase in which 

HQ 

;he discrete symmetry of the Lagrangian under ^ — is broken pj, Loinaz and Willey 
sl have used Monte Carlo simulations to calculate the critical value of the coupling constant 

that separates the two phases of the theory. 

In this work we perform similar calculations, discretizing the Euclidean quantum field 

theory (EQFT) of Eqn. [1] in terms of the two dimensionless lattice parameters 

A = Xa' fil = fiy, (2) 

where a > is the lattice spacing. (In two dimensions, both A and /ig have mass dimension 
[A] = [/Xq] = 2.) The lattice action that regularizes Eqn. [T]is 

- d 



iy=l 



(3) 



where e^, is the unit vector in the u direction. The EQFT is the continuum limit a — ^> of 
this lattice model. 

In two dimensions, the field strength and self-coupling renormalization factors and Zx 
are finite, and do not affect the phase structure of the theory. However, there is an infinite 
mass renormalization, which requires that the bare mass parameter be tuned to infinity 
as the continuum limit is taken, nl ~ /i^ln(l/a), where /i^ is the finite renormalized mass 
squared. Since A is independent of a and /Xg diverges only logarithmically as a — > 0, both A 
and /Xq vanish in the continuum limit a —>■ 0. Taking the continuum limit therefore reduces 
the number of independent dimensionless parameters from two to one, which we take to be 
the dimensionless coupling constant / = A//i^. 

We can parametrize the mass renormalization as 

nl = - Sj/^ (4) 
£^ = l(V0)^ + l/xV + ^0^-^5/xV, (5) 



2 
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FIG. 1: The only divergent Feynman diagram in (p2 theory. 

where fi^ and the finite part of S/i^ depend on the choice of renormahzation condition. 
We want to choose a renormahzation scheme in which the effective couphng constant / 
distinguishes between the two phases of the theory, which is not the case for several popular 
renormahzation conditions . We will achieve this by choosing the mass renormahzation to 
be equivalent to normal-ordering the interaction in the interaction picture in the symmetry 
phase. 

There is only one divergent Feynman diagram in 02 theory, Fig. [H which involves the 
integral 

^ N N ^ 

= ^ E E 



A^2 z_. Z_. ^2 + 4 sin2(7rA;i/iV) + 4 sm^{TTk2/N) 
d^p 1 



in the continuum limit. From Eqns. [1] and HI 



(6) 



= SXA^2 - dix^ + two-loop. (8) 

Therefore the renormahzation condition 

5i? = 3XA^2 (9) 

removes all ultraviolet divergence from the perturbation series based on the renormalized 
parametrization of Eqn. |H 

Applying this renormahzation condition. 



1(V0)^ + 1/.V + ^:0V, (10) 



dropping a constant piece in the second line. In terms of / = A//i^, the first line of Eqn. [TO] 
can be written as 



,2 



= liV^r + - 3fA,.)<p' + (11) 



On the lattice (a > 0), A^2 is finite, so we can argue that for / sufficiently small, the 
exact effective potential has a single minimum at (0) = 0. The coefficient of 0^ in Eqn. [11] 
is negative for large /, suggesting a transition to the broken symmetry phase. However, the 
effectivepotential need not be well approximated by its tree-level form at strong coupling. 
Chang [l| has shown that this transition does occur, using a duality transformation from the 
strong coupling regime of Eqn. [11] to a weakly-coupled theory normal-ordered with respect 
to the vacuum of the broken symmetry phase. 

We proceed by using Monte Carlo simulations to map the critical line in the (/Zq, A) plane. 
We determine critical values of /ioc(-^) various A, calculating the infinite-lattice-size limit 
of Monte Carlo data measured on lattices of finite size. We then impose our renormalization 
condition 

A' = Ao + 3AA^2 (12) 
using the integral representation of A^2 in the infinite- volume limit, 

POO 

A^2= dtexp[-fih]{exp{-2t)Io{2t)f . (13) 
Jo 

Here /q is a modified Bessel function of the first kind. 

For fixed A 7^ 0, we solve Eqns. [T2] and [T3] numerically to determine /i^ from /tg^. We then 
extrapolate A ^ to obtain the critical coupling constant 

in the continuum limit. We will see that this extrapolation has a nonlinear form. 



II. SIMULATIONS 



We performed Monte Carlo simulations based on the lattice action of Eqn. [3] on x 
lattices with = 32, 64, 128, 256, 512, and 1024. For each N, we set A = 1.0, 0.7, 0.5, 0.25, 
0.1, 0.05, 0.03, 0.02, and 0.01, and for each (A^, A) scanned in fi^ beginning in the symmetric 
phase and ending in the broken symmetry phase. To further constrain the data at small A, 
we performed additional simulations at A^ = 600 and 1200 for A = 0.05, 0.03, 0.02, and 0.01. 

To reduce critical slowing down, our simulations execute a Wolff cluster algorithm jj] up- 
date on the embedded Ising model after every five random sweeps of the lattice with standard 
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Metropolis updating, as in . After an initial thermalization of 2^^-2^^ Metropolis- Wolff 
cycles, we measured lattice quantities following each of an additional 2^^-2^^ cycles. 

Since these measurements are not independent, we also calculated the autocorrelation 
time T for each (A^, A, /Xq) simulation and incorporated it into our analysis. Typical autocor- 
relation times are around ten measurements, with maximum autocorrelation times around 
100 measurements for A ^ 1 on small lattices. In every simulation the thermalization time 
exceeded lOOr and we took at least 100 statistically independent measurements. As a result, 
our statistical uncertainties are quite small in comparison to systematic uncertainties. 

We use three diagnostics to determine the critical value of fi^^ where the phase transition 
occurs for fixed A. The first is the familiar peak in the susceptibility x (0^) ~ (101)^) with 
uncertainty extracted from the full width of the peak at half its maximum value (FWHM). 

The second diagnostic is the bimodality -B(/ig), a parameterization of the shape of the 
histogram of the values of measured during each simulation with fixed (A^, A, /xq) j^. 
Fig. [2] illustrates these histograms in the two phases of 02 theory: in the symmetric phase 
the histogram has a single peak around (0) = 0, while in the broken symmetry phase it has 
two peaks, around ± (|0|) 7^ 0. Constructing the histogram with an odd number of bins, we 
define the bimodality as 

B = l-^, (15) 

^max 

where hq is the number of measurements in the central bin around zero, and Umax is the 
largest number in any bin. In the symmetric phase, B <^ 1, while in the broken symmetry 
phase B ^1 (cf. Fig.E]). 

Since B depends on the specific evolution of the system, it can vary considerably for 
similar values of fi^, particularly in the symmetric phase. To smooth out this jitter, we 
consider the three-point running average B{fil) of B{fil) over /Iq, 

Bif^l) = [Bifil - Afil) + Bif^l) + B{fil + A^D] /3. (16) 

Fig. [3] illustrates the benefits of this smoothing procedure. We take as the phase transition 
point the value of /Xq for which B{(i^) is closest to 0.5, with bounds given by the fi^ most 
distant from this critical /tgc for which 0.1 < i? < 0.95 (cf. Fig. [3]). These conventions produce 
results consistent with those from the susceptibility, with comparable (though generally 
smaller) uncertainties, as shown in Table [B 
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Symmetric phase, B=0, 
(N, A, yu^)=(32, 0.05, -0.075) 



Broken phase, B=0.99, 
(N, A, fil)=02, 0.05, -0.11) 



FIG. 2: Histograms of (p for simulations with = 32 and A = 0.05, in the symmetric phase (left, 



-0.075) and broken symmetry phase (right, /Xq = —0.11 
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FIG. 3: Bimodality plotted against /ig for simulations with = 64 and A = 0.5, before (left) and 
after (right) smoothing. 

To verify that the bimodality is a robust indicator of the phase transition, we checked 
its behavior in the well-understood two-dimensional Ising model. Using the conventions 
stated above, we found that the critical fi^^ indicated by the bimodality agrees well with 
that indicated by the susceptibility in this case as well, with comparable uncertainties. Both 
observables agree with the exact analytic result. 

Finally, we extract a third estimate of the critical /Iq^ using the Binder cumulant 



u 



3(02)2- 



(17) 



For A fixed, U has a fixed point at the critical /Xq^ for any value of the lattice size N. We take 
as the critical /xq^ the value of fi^ at which U for the three largest N are closest together, 
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TARTF T- Prlflnnl r/2 


iiUlli cojCII pilcioc LidllDlLlUll lllCllCdLUl . 




A 


Susceptibility 


Bimodality 


Cumulant 


1.00 


-1.27233(16) 


-1.27258(10) 


-1.27260(45) 


0.70 


-0.95151(25) 


-0.95152(7) 


-0.95180(40) 


0.50 


-0.72080(11) 


-0.72131(9) 


-0.72130(30) 


0.25 


-0.40346(18) 


-0.40373(6) 


-0.40390(20) 


0.10 


-0.18424(11) 


-0.18432(9) 


-0.18430(20) 


0.05 


-0.10060(5) 


-0.10071(4) 


-0.10100(35) 


0.03 


-0.06410(4) 


-0.06414(5) 


-0.06420(15) 


0.02 


-0.04464(3) 


-0.04468(5) 


-0.04500(30) 


0.01 


-0.02397(6) 


-0.02399(5) 


-0.02410(10) 



with bounds given by the fi^ at which all three separate. 

This analysis of the cumulant in Eqn. [T7| produces a single critical /Xq^ for each A, while 
we have susceptibility and bimodality data for each (A^, A). Performing a linear regression 
to find the ^ oo limit of the susceptibility and bimodality data with A fixed gives us a 
total of three independent indicators of the critical fiQ^ for each A. 

We find all three values for each A consistent with each other, with comparable uncertain- 
ties (Table [T]). Combining them produces the second column in Table HTl The third column in 
Table HI] holds the corresponding critical renormalized /t^ determined from Eqns. [T^andfT^ 
while the fourth presents the values of the critical coupling A//i^ which are to be extrapolated 
to the a ^ continuum limit. 



III. ANALYSIS 



Fig. m plots the values of A//i^ in the fourth column of Table HTl and clearly rules out a 
linear A ^ extrapolation like that performed in 3]. 

Analytic investigations into the structure of scalar field theories, and super-renormalizable 
theories more generally, long ago established that correlation functions in these theories 
7Y;pically depend on logarithms of the coupling Jackiw and Templeton 



id] explicitly demonstrated the presence of such logarithmic terms in a simple 0| model. 
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TABLE II: Critical /ig^, /i^ and for different A. 



A 






A/Ac' 


1.00 


-1.27251(16) 


0.097320(46) 


10.275(5) 


0.70 


-0.95153(16) 


0.068462(45) 


10.225(7) 


0.50 


-0.72112(11) 


0.048884(32) 


10.228(7) 


0.25 


-0.40372(9) 


0.024176(26) 


10.341(11) 


0.10 


-0.18429(8) 


0.009476(23) 


10.553(26) 


0.05 


-0.10067(12) 


0.004679(33) 


10.686(76) 


0.03 


-0.06412(5) 


0.002794(15) 


10.737(59) 


0.02 


-0.04466(10) 


0.001870(28) 


10.695(163) 


0.01 


-0.02400(4) 


0.000931(12) 


10.739(138) 
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FIG. 4: Critical coupling constant A//i^ plotted against A. 

using a truncated Bethe-Salpeter equation. This non-analytic dependence on the coupling 
appears generically in more complicated super-renormalizable theories as well, including the 
02 theory we consider here. 

We can numerically investigate the effect of such logarithmic dependence by fitting the 
data in Table Hi] and Fig. H] to a function of the form 



A//}.^ = Co + ciA + C2A In A. 



(18) 



The constant Cq is exactly the continuum critical coupling constant /c we wish to determine. 
Performing this fit, we find cq = /c = 10.78(3), with = 1.21 per degree of freedom {dof). 
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TABLE III: A,/i^ ^ extrapolations of A//i^ vs. A. 



Form of A//i^ fit fc X^/dof 

fc + ciX 10.31(6) 48 

fc + ciA + C2A2 10.60(5) 5.8 

/c + ciA + C2AlnA 10.78(3) 1.2 

/c + ciA + caAlnA + csA^ 10.89(2) 0.16 

/c + ciA + C2AlnA + C3A2lnA 10.87(2) 0.13 



Performing fits with additional terms (csA^ or csA^lnA) results in even larger fc ~ 10.9 with 
very small x^/dof ~ 0.15 (Table UTTl) . Fits that do not include a term logarithmic in A are 
poor, with x^/dof ^ 1. 

We can check the consistency of these results by fitting function of A and extracting 

fc from the coefficient of the term linear in A, 

fil = d^ + \/fc + 0(\^). 

Fitting 

fil = do + X/fc + c/iP + d2\'^ In A, (19) 

we find fc = 10.77(4) with x^/dof = 1.1. As above, including additional terms in the fit 
raises fc while dramatically lowering the x^/dof, while fits without any logarithmic term 
have x^/dof > 1 (Table CV])- 

Since /t^ — > as A — 0, we should find the constant term c/q ~ in these fits, and we 
can also perform fits with d^ explicitly set to zero as an additional check. Of the fits listed 
in Table IIVI only that of Eqn. [19] has d^ vanish within uncertainty, although do is within 2a 
of zero for the fit form 

fil = do + X/fc + c?iP + daA^ 
as well. The value of fc extracted from the fit 

jj^l = X/fc + diX^ + diX^ In A, (20) 

is fc = 10.79(3) with x^ /dof = 1.0, in agreement with the values from Eqns. [T8land[T9l 

Clearly, systematic errors, particularly the choice of continuum extrapolation form, dom- 
inate over statistical errors. Tables UTTl and HVl summarize fc for various linear and nonlinear 
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TA RT TT I'M- \ 


-> extrapolations of jx^ vs. A. 




Form of /i^ fit 


fc 


xVdof 


do + X/fc 


10.24(2) 


28 


do + A//, + diX^ 


10.24(6) 


33 


do + \/fc + di\^ + d2>? 


10.55(5) 


4.1 


do + \/fc + di}? + d2>? In A 


10.77(4) 


1.1 


do + A//c + di>? + daA^ In A + d^X^ 


10.98(2) 


0.04 


do + A//c + di>? + daA^ In A + dsA^ In A 


10.93(2) 


0.05 


A//c 


10.27(2) 


47 


A//c + diA2 


10.31(7) 


49 


A//c + «iA + a2\ 


10.ol(oj 


5.5 


\/fc + di\^+d2>? In A 


10.79(3) 


1.0 


A//c + diA2 + d2A2 In A + diX^ 


10.90(2) 


0.18 


X/fc + di}? + d2>? InA + dsA^ In A 


10.89(2) 


0.14 



extrapolations, along with the goodness of the fits, /dof . Neglecting fits with /dof ^ 1, 
we adopt a final result of 

fc = lO.Stf, (21) 
to be consistent with the numbers in Tables IIIII and IIVI 



IV. DISCUSSION 

Since our approach closely parallels that of Loinaz and Willey j^, it is distressing that 
our final result disagrees so strongly with the fc = 10.261q4 reported there. However, our in- 
dividual data points are largely consistent with theirs, as shown in Fig. [51 The disagreement 
between our final results comes almost entirely from the nonlinear continuum extrapolations 
discussed above. 

Re-analyzing the data in , we find (Tables |V] and IVip that they are consistent with 
all the nonlinear fits considered in our analysis above. Both linear and nonlinear fits all 
have x^/dof ~ 0.5. While nonlinear extrapolations were not required by the data in [3], by 
considering only the linear case Loinaz and Willey overlooked significant systematic effects 
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FIG. 5: Our data for A//i^ compared with the results of 3| (empty circles). 



TABLE V: A,/i^ extrapolations of X/fii vs. A, for 



a. 



Form of X/fil fit 



fc 



xVdof 



fc + ciA 

fc + CiA + C2\^ 

fc + ciA + C2AlnA 
fc + ciA + C2AlnA + csA^ 
/c + ci A + C2A In A + csA^ In A 



10.32(5) 
10.34(9) 
10.34(12) 
10.18(25) 
10.23(22) 



0.39 
0.49 
0.51 
0.59 
0.63 



due to fit form. In Tables |V] and IVH we see 10.2 < fc < 11-9 in fits with x^/dof < 1, 

.10 

.05- 



consistent with our result fc = 10. 8^'^° 



Several other authors have also calculated the critical coupling fc in 02 theory using a 
variety of methods, numerical schemes, and analytic approximations. These approaches, 
summarized in Table IViTt produce a large spread of results, of which ours is the largest. 

The density matrix renormalization group result fc = 9.982(2) |l^ is notable for its ex- 
tremely small claimed uncertainty. However, this result relies on linear A ^ extrapolations 
like those in [sl, performed with just two data points at A = 0.1, 0.25, for fixed lattice size 
N = 500 or 1000. A linear 1/A^ extrapolation is then performed using the two resulting 
values. Thus, we expect this result to suffer from the difficulties discussed above. 

The diffusion Monte Carlo result fc = 10 ± 0.8 ± 0.4 [15] agrees with our result within 
its relatively large statistical and systematic uncertainties. The Gaussian effective potential 
results fc = 10.272 12|, fc = 10.211 (J, and fc = 10.21 Q (the last of which coincides 
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TABLE VI: A — > extrapolations of /i^ vs. A, for 



-2 

Form of /i^ fit 


Jc 




"0 + Jc 


iU. Zo^o j 


U.DO 




1 n 'i7(0i\ 
iu.o < i^yj 


n AA 




lU.fJ 1 y J-U ) 


n "^9 






n 9"? 


do + A/ Jc + fllA + Ci2A in A + (13 A 




n 09 

U.UZ 


UQ -\- Aj Jc ~r Oj\A + U2^ inA + UsA m A 


11 fi/l /'I 0"\ 


n ni 


\ / f 


in 9/1 ('9^ 


n ^1 
u.o / 


A/Zc + diA^ 


10.32(5) 


0.38 


A//c + diA^ + d2A3 


10.35(9) 


0.49 


A//c + diA2 + d2A2 In A 


10.35(12) 


0.50 


A//c + diA2 + d2A2 In A + diX^ 


10.19(25) 


0.59 


\/fc + di}? + d2>? InA + dsA^ In A 


10.24(22) 


0.63 



with the oscillator representation result) are the next closest. Both the Gaussian effective 
potential and oscillator representation methods incorrectly predict a first-order phase tran- 
sition, in violation of the Simon- Griffiths theorem 2^, which requires the 02 theory phase 
transition to be second order. 



V. CONCLUSIONS 

We have used Monte Carlo simulations to obtain an accurate lattice measurement of the 
continuum critical coupling constant fc = 10.S~^'q2^ in 02 theory, improving the previously 
reported Monte Carlo result yl. 

n 

While our data are consistent with those reported in [3[, our improved precision forces 
a nonlinear extrapolation to the continuum limit, producing a significantly different contin- 
uum result. The older data are compatible with these nonlinear extrapolations, although 
such nonlinearity was neither required nor investigated previously. Applying nonlinear ex- 
trapolations to the older data, we obtain continuum results consistent with our own. 

Significantly, nonlinearity - in particular terms logarithmic in the lattice coupling A - 
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TABLE VII: Various results for the critical coupling 



Method 


Result 


Reference 


Monte Carlo 




This work 


Gaussian effective potential 


10.272 


[12j 


Gaussian effective potential 


10.211 


yy 


GEP and oscillator rep. 


10.21 


[13j 


Spherical field theory 


10.05 


M 


Diffusion Monte Carlo 


10 ±0.8 ±0.4 


[15] 


Density matrix RG 


9.982(2) 


[16j 


Continuum light-front 


9.91 


m 


Connected Green function 


9.784 


[12j 


Coupled cluster expansion 


3.80 < fc< 8.60 


[18j 


Discretized light-front 


7.325, 7.71 


[19] 


Discretized light-front 


7.316, 5.500 


[20, 21J 


Random phase approximation 


7.2 


[22j 


INTon-Gjiimman AAanatinn^il 


6.88 


[23] 


is expected analytically. This convergence of 


analytic theory and numerical dat 


a provides 


additional evidence that our improved result 


is accurate and reliable, and can be used to 


evaluate analytic approximations. 
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